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Abstract 

We analyze the singlet axial form factor of the proton for small momen- 
tum transferred in the framework of QCD sum rules using the interpolating 
nucleon current which explicitly accounts for the gluonic degrees of free- 
dom. As the result we come to the quantitative description of the singlet 
axial constant. 



1. Introduction. 

The investigation of polarized deep inelastic scattering is one of the 
most attractive field for theoretical consideration since it provides an im- 
portant insight into the structure of hadrons and opens a large area of 
subtle dynamical phenomena associated with the spin dependent case. In 
the last several years there has been an increasing interest in the deep in- 
elastic structure function gi{x). It was provoked by the EMC result on 
the scattering of the longitudinally polarized muon beam on a longitudi- 
nally polarized hadron target. The unexpectedly small asymmetry found 
by EMC has led to the so called "spin crisis in the parton model" and has 
raised a number of questions of understanding the dynamics of the proton 
spin on the parton level, namely, how the nucleon spin is build up from the 
spins of its constituents. An enormous flood of theoretical investigations 
was generated in order to resolve the current "spin problem" 

The EMC measurement of the first moment of the polarized structure 
function T\ can be interpreted, via the Ellis- Jaffe sum rule 

T\{Q') = j'^dxgl{x, Q') 

4{(a!f(0).i5C-(0))(l-(^)-3.5833(^)^20.:53(^)^) 

.5a-(0,Q^)(l-l(^)-1.0959(^)^)}. (1) 

as a first ever measurement of the singlet axial constant G^^(O) of the 
proton. The last one turns out to be unexpectedly small in contradiction 
with the naive parton model where it is fairly close to unity. The EMC 
reported the result for G'^^(O) which is compatible with zero. The new 
experiments are performed to check their measurement of gi{x) and to 
measure an analogous neutron function gi{x). The recent analysis of 
proton and deuteron data gives G^^\o) varying from 0.20 ±0.11 to 0.36 ± 
0.05, that is still far from unity. So the problem reduces to the evaluation 
of G^^(O) because the other two axial constants can be extracted reliably 
from the data on deuteron (or He^) target and on neutron and hyperon 
/3-decays. In this paper we calculate it in the framework of QCD sum 
rule approach which till now seems to be the most powerful method for 
extraction of information about the low energy properties of hadrons and 
the closest one to the first principles of the theory. 
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In eq.(|I|) the functions G^XiQ'^) 

are form factors at zero momentum 
transferred in the proton matrix elements of axial currents 

{N[p2A2)\jllmN{piAi)) 

= 4^\p2) {G^hQ'h,l6 - Gf{Q')q,^,) (2) 

where ? is a SU{3)f index, q = P2 — Pi and = —q^. There is an 
important difference in the behaviour of induced pseudoscalar form factors 
at small momentum q. Here, the singlet pseudoscalar form factor does not 
acquire a Goldstone pole at = 0, even in the chiral limit, contrary to 
the matrix elements of the octet currents. It is known that this limit, in 
which the masses of the three light quark flavours are neglected, is not 
far away from the real world of hadrons. In this limit, there exist eight 
massless pseudoscalar mesons serving as Goldstone bosons. However, the 
ninth pseudoscalar, the ?7'-meson, remains massive. In the following this 
property will be used to extract a value of G^^\o) from the sum rules. 

It has been established ^ that the first moment does not measure 
the contribution of the quark spins to the proton one. This happens due 
to the anomalous nonconservation of the singlet axial current. For this 
reason, we display the role of this profound feature of the theory from the 
very beginning exploiting the equation for the anomalous divergencel^: 

d,ji'J = ^ii: m,qi,q " (^) G%G%, (3) 

where Nf is a number of flavours (later Nf = 3). Taking the divergence of 
eq.(|2|) for the singlet axial current and making use of the last expression 
we come to the relation which directly connects, in the chiral, limit the 
nonforward matrix elements of the gluon operator to the effective form 
factor 2mNG^^]j{Q'^) = 2mNGf{Q'^) + Q^Gf{Q'^) that is equal to the 

2mNGf{0) at = 0. 



^Throughout the paper, we adopt the conventions in Itzykson and Zuber 



2. Effective axial form factor in QCD sum rules. 

For a long time all calculations of the nucleon characteristics use a par- 
ticular three-quark current introduced by loffe [)5|. When one makes an 
attempt to evaluate the matrix elements of quark-gluon or gluon opera- 
tors, one faces evident calculational difficulties, moreover, the final sum 
rules are aggravated by extra UV logarithms due to mixing of operators 
and, therefore, the calculations are affected by noncontrollable uncertain- 
ties 

In field theory, the usual statement that the nucleon mainly consists of 
three quarks means that the 3 quarks 3 quarks Green function (three 
quarks are in a state with nucleon quantum numbers) has a pole at the 
mass of the nucleon, with the total angular momentum J = |, with a 
significant residue. The fact that the nucleon is not just a three quark 
state means that the nucleon pole also occurs, albeit with smaller residue, 
in Green functions such as 3 quarks+g 3 quarks+g. For this reason, one 
is forced to introduce a more sophisticated interpolating proton field which 
explicitly contains the gluonic degrees of freedom: 

tjg{x) = e'^^ {u\x)C^,u^{x)) ^,j,a,p{gGlp{x)t-d{x))\ (4) 

The latter was investigated in ref. [j7| and checked in the calculation of pro- 
ton gluonic form factor normalized to the fraction of nucleon momentum 
carried by gluons. Recently, making use of this current the twist-3 and 
twist-4 corrections to Bjorken and Ellis- Jaffe sum rule have been found 
0. The advantages of this current are straightforward: the calculations 
are drastically simplified, sum rules are free from additional divergences 
that are not removed by the single Borel transformation. At the same 
time, the applicability of non-dimensional regularization may be spoiled 
by the power UV divergencies, appearing due to the high mass dimnen- 
sion of this current. For the same reason, sum rules may be affected by 
the vacuum condensates of higher dimensions, reducing their outcome to 
semiquantitative estimates. 

The usual technique of QCD sum rules is to extract the nucleon ma- 
trix element of local operator from the appropriate three-point correlation 
function. This correlator is the sum of different tensor structures each 
characterized by the relevant invariant amplitude W^'^\pl,p2.,q'^). 
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(5) 

In practical calculation it is advantageous to consider W^^\pl,pl, q^) (here- 
after referred to as W) because of its lower dimensionality. Another reason 
in favour of this choice is that it does not lead to the fictitious kinemat- 
ical singularities in q^ as the last term in eq.(^ does. For this invariant 
amplitude we can write the double dispersion representation 

^{Pl:P2^q ) = — L L dSidS2- 2V + ' ' ' ' ^ 

TT^ JO Jo - pf^(^s2 - pi) 

where the ellipses stand for the polynomials in pi and P2 which die out after 
the double Borel transformation has been applied. For the physical spectral 
density we have accepted the conventional "resonance plus continuum" 
model: 

^pcont{sij S2, Q^)(l - 9{ao - si)9{ao - S2)). (7) 

The function in front of the double-pole term is a combination of form 
factors we are interested in up to certain overlap Xq between the state 
created from the vacuum by ijg and the nucleon state 

{0\r]GmNipA))=m%XGU^N\p)- (8) 



So, our aim is the evaluation of the correlation function (|^) in QCD. In the 
case when all the momenta {—Pi) ~ {—pD sufficiently large (of an 

order of IGeV^), the leading contribution comes from the domain where 
all distances are small. Thus, the standard machinery of short distance 
expansion are applicable, allowing one to express the final result in terms 
of quark and gluon condensates. The problem is modified drastically if 
the squared momentum transferred becomes small {Q^ <C {—pf)) because 
the relevant t-channel distances can be large. In this case the OPE has a 
twofold structure 0. Terms of the first type arise from the SD(I)-region 
when all intervals ~ ~ {x — are small. Another contribution 
comes from SD(II)-region (bilocal power correction) which originates from 
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the distances ~ (x — yY- The necessity for the bilocal power 

corrections can be traced from the fact that the ordinary QCD Feynman 
diagrams contributing to the form factor at moderately large in the 
limit of small possess logarithmic non-analyticities {Q'^)^lnQ'^, which 
signals that large distances come into play [T^]. Therefore, we have to 
subtract such a perturbative behaviour from the corresponding graphs and 
add the "exact" correlators which account for the nonperturbative effects 
and thus possess the correct analytical properties as goes to zero. So 
the OPE in the case when the momentum transferred can be arbitrary 
small has a modified form [TT| 

= T.CfL(,)(p\,pW){Oi) + Y.j d'xe""C%^„^(x^)Wr(:c,ql (9) 

d i 

where, as was mentioned above, the coefficients C^sD{i)iPiiP2^ regular 
in the limit ^ 0. The second term determines the large t-distance 
contribution. Here W^^ are the two-point correlators 

W^\x, q) = J d'ye^^^{0\T {g{y)0,{x, 0)} |0) (10) 
of operator in question Q{y) = [^^^G'tj^{y)G^ {y) and some nonlocal string 



47r ^ll 

operator with definite twist (not dimension) [jT^, [TBj that arises from the 
OPE of the T-product of nucleon currents: 

T{m{x)nGm = EC^E(//)(^')0.(^,o) (11) 

i 

The bilocal power corrections cannot be directly calculated in perturbation 
theory but we can write down the dispersion relation for them 

Wr(:c,q)='-rd /'^''^^Y^ , (12) 

TT "'0 s — 

assuming the standard spectral density model with continuum to start at 
some threshold sq and finding in some way its parameters. We always do 
this constructing auxiliary sum rules. There is no need in additional sub- 
tractions in eq.(|T2|) because one always deals with the difference between 
the "exact" bilocal and its perturbative part; so due to the coincidence of 
their UV behaviours the subtraction terms cancel in this difference. 

To simplify the calculation of the local power corrections, it is convenient 
to use fixed-point gauge for the background gluon field {x — xo)^B"^{x) = 0. 
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We chose the fixed point in the vertex of the gluon operator xq = 0. The 
quark and gluon propagators in this gauge up to the order 0{G) looks like 

m 



^jyv^5 y) 

.r(f-l) g^J-b 



. d 

47r2 



V(-- 



167r2 



■A2] 



87ri 



[-A2]^ 



S{x,y) 

r(|) ^ 



27rf 



-A2] 



r(|- 



A, 



-A^ 



Til) ^ 



[-A2] 



d 1 



(13) 



where A = x - y, G»t = ff/^'G^, „ . . „ 

gt°'G^j^ for the quark one, the generators are normalized by Sp{t°'t^ 
For the noncoUinear quark condensate we use the following expansion in 
terms of local vacuum expectation values 



y, = (//"^^G^j, for the gluon propagator and G^^y = 

2" • 



1 



43 



ml{il)il)) 



l-(J 



(3a 



(14) 



The ellipses stand for the higher dimension vacuum condensates. 

Calculating the diagrams depicted in the first row of fig. 1 we come to 
the Borel sum rule for effective axial form factor at moderate values of the 
momentum transferred with the Borel parameters ri and T2'. 



T1T2 



TT^ JO JO TtTo 
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18 L [Ti+T2y [ri+T2y 

^ ^2L_(lll2)^r I r>2 (^1^2) J 1 

144 L (ri+r2)'^ (ri+r2)^ J 
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32 [Ti+r2y 



T I r.2 '^^^^2) J 
ri+T2)'^ [ri+T2)^ 



J ^ 2 

J12 - -mo 



(rir2] 



Ti+r2)' 
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(15) 



6 



where the continuum double spectral density is 

5 

Pcont{si, S2, Q^) = E S2, Q^), (16) 

and each term in a sum is found from the corresponding diagram in fig.l 
(appropriate formulae are given in appendix B): 

.a)(....g^) = S(f)^..)g^(i-^)^(2 + ^), (17) 

P(2)(.i,..,Q^) = §(f)-?MQ^{^(i-^)\2- 



2\ Rij ^ rI ^' 



(18) 



/3(3)(S1,S2,Q^) 

(19) 

g^) = -'^(^f'"li--)Qf - jrf (2 + ^) . (20) 

/3(5)(S1,S2,Q^) 

^-f(^)^«.^MQ^{^(i-|.)t2.|.)4(i-|.).g^^), 

(21) 

and 

(J = Sl + S2 + Q^ i?(si,S2,Q^) = (J^ - 4siS2. (22) 

The functions are originated from the diagrams in the first row in 
fig.l and are given by the following expression: 

Jnm{r^, Q') = /' dxx^-'x^exp { ---^] (23) 



We state that contrary to the refs.[§j where the sum rules with the 
same interpolating nucleon field were dominated by the contribution from 



the highest dimension operators, our sum rule is not affected by them: 
the coefficient functions that are determined to the leading accuracy by 
tree and one-loop diagrams vanish identically. Therefore, we do not meet 
the problem of breakdown of OPE for the correlator in question. The 
absence of higher condensates contribution unsuppressed by a number of 
loops is directly connected with chiral structure of the interpolating nucleon 
field and the tensor structure chosen for investigation in the three-point 
correlation function. 

As was mentioned previously one could not put = in the eq. 
(fTBl) because though finite it is including contributions non-analytic in this 
point. This is typical example of the mass singularities. Therefore, follow- 
ing the method outlined in the middle of this section one should subtract 
the perturbative behaviour from corresponding graphs (the diagrammatic 
representation for them is shown in the second row of fig. 1, while the 
explicit expressions are written in the appendix A) and add the terms with 
correct singular structure in . It is clear that the singularity should be 
located at the threshold of the first prominent resonance in the correspond- 
ing channel. 

3. Bilocal corrections. 

The simplest bilocal correction (first picture on fig. 2) is given by the 
convolution of the coefficient function involving the quark condensate with 
the two-point correlation function of operator in question and some point- 
splitted gluon operator coming from the OPE of nucleon fields 

2 

= i (^) /'i^e'™{0|T{G°„{y)G°„(j/)(G»,(:E)G°„(0)-G°„{x)G°,{0))}|0) 

(24) 

As was mentioned earlier we cannot calculate it in perturbation the- 
ory but we can reconstruct it from the information about its large-Q^ 
behaviour. To this end we write down the dispersion relation for it of the 
type represented by eq. ([T2|) and use the standard "resonance plus con- 
tinuum" spectral density model, with Ty'-meson. It is likely to be the only 
prominent singlet pseudoscalar both in quark and gluon channels. 
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p(s) = TTmlJ^\if^^^(t)^^){xq) + f^^\xq)(t)^^}{xq)]5{s - m^,) + e{s - sq)p^^{s). 

(25) 

Here we have parametrized the matrix elements of some gluon operators 
between the vacuum and r/'-meson state as foUows: 



Nja 
47r 



= {xi^qx-xxqK)[ifn^U^ry\^q)^ f!y\^(l)(I^S\^(l)]- (26) 



(i) 

In the last line the wave functions 0y can be related in the standard 



way to the usual ones {(f^^'), describing the light-cone momentum fraction 
distribution of gluon inside meson. 



In eq. (|2^) we have kept only the leading twist wave functions which 
reproduce the leading nonanalyticity in the corresponding contribution of 
the local {uu) power correction. As will be shown below, we account for 
other contribution using simple recipe which results from our consideration 
(originally it was proposed in Ref. [|TB[). We can find the unknown overlaps 
//^*^ constructing the auxiliary sum rules. It turns out, that due to the 
antisymmetrical tensor structure involved the contribution of ordinary local 
power corrections with gluon condensates, unsuppressed by a loop factor, 
vanish identically in the theoretical part of the sum rules. For this reason 
we account for nonperturbative effects introducing the concept of nonlocal 
gluon condensate [jT?! which corresponds to infinite series of local ones. It 
can be appropriately decomposed into two tensor structures multiplied by 
corresponding form factors [1T5|, [T5j : 

{0\G-'Jx)E'\x,0)Gl{0)\0) 



{G 



2\ 



12 



dDA{x'^)^ 

~^ yg ^Lv-^ a ~i~ gpa-^pL'^v g pa-^v-^ p gvp-^ p-^cj) j ■ 1 / 

This form explicitly separates out the term proportional to Dj^a{x'^) which 
violate the abelian Bianchi identity, while the second term satisfies it. It 
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was shown that hnear confinement occurs when D]\[a{x'^) is present in (^) 
while the second term does not contribute to the string tension 



In the calculation of fjy ^ and fjy ^ constants only abehan form factor 
contribute. We present it in the form of a-representation [p^] : 



Da(^') = l"daf^{a,Xl)e-'r, (29) 
and use a (5-shaped ansatz for the distribution function A^): 

/^(«,Ai) = %-Ai), (30) 

where 1/Xa is an abelian correlation length of the vacuum fiuctuations, it 
can be expressed in terms of vacuum condensates X\ = ^g^ [uuf / [G"^) « 
O.OSGeF^ at IGeV"^ []23]. One comment concerning eq. ( |29| ) is that in 




deriving a QCD sum rule one can always perform a Wick rotation xq ixq 
and treat all the coordinates as Euclidean, < 0. 

Proceeding in the standard way we obtain the following sum rule: 

-^^is^^^^^Hii^j + T^V^^ii^l^-M^''' ^^^^ 

where 

2 

E2(x) = l-(l + x+|-)e-^ (32) 

We stress that due to the fact that one of the gluon currents has nonzero 
Lorentz spin leads to the absence of direct instantons to the polarization 
operator of interest []2T| . This property may be considered as a counterpart 
of the absence of the topological ghost pole in the formfactor of interest 
while this pole does appear in the formfactor of the conserved quark-gluon 
current, related to the low-energy nucleon structure. 

We take now the limit oo and obtain the local duality relation: 

^2 r Ai)_Njas( 1 Nja^sl Nf a^^ A 
^v'UU -^ U^^"^ + ~^~^^^^l' ^^^^ 



The value of continuum threshold is found from the requirement of the 
most stable sum rule (|3T|). Straightforward analysis gives us the value 



So = 2.5^6^^ which coincides with the one obtained in ref. 

Keeping the contribution due to the nonlocal gluon condensate would 
exceed the accuracy we are pretending to because, as it was mentioned 
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above, we do not calculate the corresponding term in the OPE with local 
power corrections which are obviously small. The leading non-zero con- 
tribution coming from non-local condensate was required to analyze the 
stability of the sum rule and to determine the continuum threshold. How- 
ever, the relative numerical value of the non-local condensate contribution 
is small (like that of the dropped local term), and we neglect it in what 
follows. 

For the overlap factor an analogous relation taken to the same 
accuracy is: 

The net form for the additional term for the axial form factor at small 
momentum transferred looks like 



3 ' ' I (47r)2 V TT 



N; g2 + 

{uu)M^e' 



Q2 j 2 



2 
3 




In the last line we have substitute the local duality relations for residue 
factors and take the asymptotical form for the wave functions 

(^(2)(c^) = A2{){a-a)a^a^ (36) 

Note that in eq. ( |35|) we have put ti = r2 = 2M'^ in order not to 
introduce the asymmetry between the initial and final states and to make 
contact with Borel parameter of the two-point nucleon sum rules. We 
can observe that at large the bilocals vanish faster than the term it 
is correcting ~ Q^. Note, that although the local power correction may 
vanish for ^ the modified version of nonanalyticities Q'^^ In alive 
in this limit. 

As can be easily seen, the leading nonanalyticity in the (ww)-term in 
the correction found and in the expression for form factor cancels. It is 
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replaced by the combination sq + which is "safe" in the limit 0. 
In the same way we may correct the other nonanalyticities. In large 
limit where the original OPE must be valid the bilocal corrections must 
be absent. As we have seen, the residues of the physical spectrum can 
be found from the requirement that the bilocal power corrections should 
vanish faster at large than the contribution they are correcting. Then, 
in general the correction term is given by the following equation (we omit 
all unnecessary constants) 



Ti/res+co7it 
BL 



dss 



n— 1 



So 

S + Q2' 



(37) 



Using this simple recipe one can easily modify all perturbative non- 
analyticities of the effective axial form factor in the small momentum trans- 
ferred limit. After all of them have been corrected properly we can take the 
limit ^ and obtain the sum rule for singlet axial constant directly. 
As was noted in Ref. []7|, one should not try to fix the parameter of con- 
tinuum threshold ctq from the sum rule with a new current. Therefore we 
take the continuum threshold usual for the sum rules involving the nucleon 
and consider it in the local duality limit. Combining all contribution we 
come to the following equation 



5 A 2 ^(0) _ 



223 



751 2d 



2732 




+ 



1 31 2 

-i?3 + ^moi?2 



1479 2 r. 

R5 + —TT^TTlnR 



21033 



0^4 



where 



Rn = 



n-1 



nm. 



n 



1 



(38) 



(39) 



and its origin was clarified above. The sum rule imply that we model the 
continuum by effective spectral density that includes all ones which are 
nonzero for s > 0. 

We use the standard ITEP values of condensates rescaled to the normal- 
ization point fi^ ~ m% ~ IGeV^ with the appropriate anomalous dimen- 
sions: a = -{2TTf{uu) = Q.QlGeV^, ml = {ug{aG)u)/ {uu) = 0.65GeV^; 
also, we use the overlap value = 2(27r)'*A^ = O.SGeV^ and continuum 
threshold ctq = 2.5GeV'^, providing the better accuracy of the calculation 
of partonic densities The value of the strong coupling constant at 

IGeV'^ is taken to be ag = 0.37 that corresponds to A = IbOMeV. We 
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obtain the following numerical value of the singlet axial constant 

^^^^(O) = 0.2. (40) 

Varying the parameters in the reasonable range will result in the variation 
of the quantity within the 50%. The main uncertainties come from the 
errors in estimation of the ^-channel continuum threshold sq and the overlap 
Xg of the nucleon state with that created by the new current. 

Due to the anomalous non-conservation of the singlet axial current the 
singlet axial constant is not a renormalization group invariant. Therefore 
in order to compare our prediction with the experimentally mensurable 
quantity, we have to evolve it from QCD sum rule scale /i ~ IGeV up 
to the one of EMC-SMC experiment which is = lOGeV^ exploiting the 
one-loop solution of RG equation: 

Gf{0, Q') = Gf{0,f,')exp!^£^^[a,{Q') - ^.(/i^)]} , (41) 

where the anomalous dimension 72 = 16Nf and as usual (3q = 11 — |A^/. 
However, the sensitivity to the QCD radiative corrections is poor until very 
large is attained and account for them would exceed the accuracy of 
our estimate. Nevertheless, the value obtained are in reasonable agreement 
with the new world average value for the singlet axial constant. 



4. Summary. 



In summary, we have calculated the singlet axial constant in the QCD 
sum rule framework for the form factor type problem at small momentum 
transferred and find the value in good correspondence with experimental 
one. We should mention that in our letter ^ G^^^ was somewhat over- 



estimated because the next-to-leading twist bilocal power corrections were 
not accounted for and the continuum was not properly subtracted. 

In ref.[)2^ the pioneering attempt was undertaken to evaluate G^^\o) 
by QCD sum rules in a way similar to the calculation of the octet axial 
constant [j25j. Due to the presence of the gluon anomaly in the induced 



vacuum condensates the problem differs significantly from the one for the 
G^^(O). This feature was incorporated in the calculation but nevertheless 
the authors did not come to the reasonable quantitative prediction of the 
singlet axial constant. It was conjectured that the OPE breaks down for the 
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singlet axial current in the axial-nucleon-nucleon vertex. At the same, we 
did not observe any evidence of the divergence of the OPE in the correlator 
under investigation: the contribution of the highest dimension vacuum 
condensates unsuppressed by a number of loops is absent. From the other 
side, the small values of the power corrections result in the good accuracy 
of the local duality approach and in the strong dependence of our result 
on the continuum threshold. 

A possible hne of development would be to estimate twist-three gluon 
contribution into the moments of the transverse spin structure function g2 
PE| and the x-dependence of the twist-two polarized gluon distribution in 
the nucleon. However, the latter would require an elaboration of the new 
procedure for separation of the large and small distances in the effective 
four-point correlator. 
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Appendix A. 



In this appendix we present some useful integrals needed in the addi- 
tional factorization of large and small distances in the ordinary Feynman 
diagrams. Keeping the lowest twists contributions into the perturbative 
parts of the bilocal correlators we find: 

.A {x{k + q)Y ^ i{-iy 1 r r(r-f) 

(42) 
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(43) 



Ji27rY [k^-lY (47r)f 2r(r)l ^-^-^ 



(47r)i 2r(r) r-f-i 



(44) 



where q = ao, and CJJ^ = —rr^ — r, are binomial coefficients. 

^ ^' n m\(n—my. 

Using these results we obtain for the factorized diagrams in the second 
row in fig.l: 



TT^ 18 V TT / [Tl+T2y ^ 



+[4(rir2) - 6(ri+r2)']L4 + ^[4(rir2) - lbin+T2f]Le)e(S2) ^ (46) 



^^^^ = i^f^l -^(^-)7M.([4(nr2) - (n+r2)^]L2 



7r2 144 V Try " {ri+Ti)^ - 
+[4(rir2) - 2(ri+r2)']L4 + ^[4(rir2) - 3(ri+r2)']L6)e(^, (47) 
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TT^ 8 V TT 



+[2(rir2) + 2(ri+r2)']L4 + -[2(rir2) + 3in+T2)'']Le)ei^ , (49) 



and 



-f'n — 



n/2 



rir2 
>2 (ri+r2)2 



)-5i(n/2-l) 



with 



(50) 



(51) 



where we have used the M^-scheme. 

As can be easily seen the logarithmic terms in eqs. (y4.4 — AS) repro- 
duce the leading non-analiticities in the expressions for ordinary diagrams 
contributing to the form factor at the moderately large Q^. Therefore this 
perturbative " long-distance" behaviour cancels exactly in the difference of 
the diagrams in the first and second rows in fig.l. 



Appendix B. 

Here we present the formulae which enable to represent the integrals 
appeared in the calculations in the form of double spectral representation. 
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Figure captions. 



Fig.l. Contribution to the effective axial form factor in the QCD sum 
rules approach. The difference of the first and second rows defines the 
SD(I)-regime. 



Fig. 2. Generic form of the bilocal power corrections entering the OPE 
with different coefficient functions. 



